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The pinning of the Fermi level to the van Hove singular- 
ity and the formation of flat bands in the two-dimensional 
t-t' Hubbard model is investigated by the renormalization 
group technique. The "van Hove" scenario of non- Fermi- 
liquid behavior for high-T c compounds can take place in a 
broad enough range of the hole concentrations. The results 
are in qualitative agreement with the recent ARPES data on 
La2Cu04- 



The nature of normal state and the mechanisms of su- 
perconductivity in copper-oxide high-T c compounds are 
still a subject of hot discussions |]-§]. A number of 
experimental data on transport properties and ARPES 
give an evidence of a non-Fermi-liquid (NFL) character 
of the normal state in underdoped regime and, as a con- 
sequence, of a strict confinement of current carriers in the 
CuC>2 planes. The latter results in an anomalously weak 
sensitivity of normal-phase and superconducting proper- 
ties to the impurities outside the planes, hoppingchar- 
acter of the transport along c-axis etc. Anderson (l| has 
put forward the idea that a 2D system can demonstrate 
a NFL behavior at arbitrarily small interelectron repul- 
sion U owing to a finite phase shift at the Fermi energy. 
In this scenario the confinement of current carriers is ex- 
plained by "quantum protection" owing to charge-spin 
separation and incoherent (non-quasiparticle) character 
of electron (or hole) motion. However, theoretical un- 
derstanding of NFL state in quasi-2D electron systems 
turned out to be very difficult. Modern field-theoretical 
investigations show that such a state, probably, does 
not occur in general 2D and quasi-2D cases. 

At the same time, peculiarities of electron spectrum 
can lead to NFL behavior. In the presence of van Hove 
singularities (VHS) near the Fermi level the marginal 
Fermi-liquid B or NFL behavior |)| can be naturally 
derived. Already in the leading order of the pertur- 
bation theory in U the marginal dependences of the 
electron scattering rate and specific heat take place, 
Im£(kF,e) oc |e|ln|l/e| at |e| 3> \fi\ and SC cx 
Tln 3 (t/ max(|/i|, T)), and the resistivity demonstrates 



the behavior p oc Tln 2 (i/max(|^|,T)) (lOj (T is the tem- 
perature, e is the energy, fi is the chemical potential cal- 
culated from VHS, and t is the hopping integral). Close- 
ness of the VHS to the Fermi level results in a ID-like 
behavior of the perturbation expansion pd| , but does not 
lead to the Luttinger-liquid fixed point since the system 
turns out to be unstable with respect to formation of 
magnetic or superconducting ground state. It has been 



shown []l l| |l 3|] in the framework of 2D t-t' Hubbard model 
that the closeness of the Fermi surface to VHS leads to 
the instabilities of the Fermi-liquid state in the finite 
range of the concentrations near the van Hove filling. 
The type of the instability (magnetic or superconduct- 
ing) depends on the value of t'/t. 

This "van Hove scenario" (for a review of early consid- 
erations, see, e.g., Ref. Q) seems to be very attractive 
since various data demonstrate the closeness of VHS to 
the Fermi level in high-T c compounds at optimal dop- 
ing or pressure. Recent ARPES data demonstrate that 
the Fermi level of La2- a: Sr a ;Cu04 is close to VHS in the 
range x = 0.20 -j- 0.30 |^5|, which corresponds to over- 
doped regime. One of the two bands, that occur due to 
bilayer splitting (antibonding band) in Bi2212 system, is 
also close to VHS 16 . A significant and still unsolved 



question is why the Fermi level should be close to VHS 
in a broad doping range. Although some considerations 
were performed (l^^^], they did not take into account 
the renormalization of electron spectrum, e.g., the for- 
mation of the flat part of the spectrum, which should 
change substantially the results of these approaches. 

The flattening of the spectrum was first studied theo- 
retically for a single hole in an antifcrromagnetic back- 
ground Hi. The flat region leads to some peculiarities 
of electron properties near the metal-insulator transition 
flisfl , in particular to the pinning of the chemical poten- 
tial. While above-mentioned approaches describe cor- 
rectly the situation at small hole doping (the doped Mott 
insulator regime) , these approaches meet with difficulties 
near the optimal doping where the concentration of holes 
is not small. In this Letter we argue that besides anti- 
ferromagnetic fluctuations there is another factor which 
results in occurrence of the flat part of electron spec- 
trum, namely, the critical fluctuations that occur near 
the van Hove band filling. We also reconsider the issue 
of the pinning of Fermi surface to van Hove singulari- 
ties by proposing a new scheme of renormalizations of 
the energy- and momentum-dependent electron Green's 
function. We show that the pinning of the Fermi surface 
to VHS is an universal feature of 2D systems which is 
connected with the formation of a flat region of electron 
spectrum near the (tt, 0) point. Both the phenomena, the 
pinning and flattening, will be described within the same 
renormalization group (RG) approach. 

We start from the t-t' Hubbard model on the square 
lattice, 



(1) 



with the bare electron spectrum 



1 



£k = — 2i(cos k x + cos k y ) — 4t'(cos k x cos k y + 1) — fi 

(2) 

Hereafter we assume t > 0, t' < (which is the case 
for hole-doped systems), < \t'\/t < 1/2. For arbitrary 
t'/t the spectrum (Q) contains VH singularities connected 
with the points A = (7r, 0), B = (0,7r). The chemical 
potential /i measures the distance between VHS and the 
Fermi energy, so that at fi = VHS lie at the Fermi level. 
We suppose that t'/t is not too small, and therefore the 
Fermi surface is not nested. It was argued ||[l2|,|l3| that 
provided that the band filling is close to van Hove one, the 
vicinity of the van Hove points gives a dominant (most 
singular) contribution to electronic and magnetic prop- 
erties. Thus we can use the expansion of the spectrum 
near A(B) van Hove points 



2i(sin 2 ipk x — cos 2 <pk 2 ) — p, (3a) 
sin 2 (fiky) — p (3b) 



el = 2t(cos 2 ipkl 



where k x = tt — k x , k y — tt — k y , ip is the half of the 
angle between asymptotes at VHS, 2ip = cos _1 (— 2t'/t). 
The chemical potential p is obtained from the electron 
concentration n — J . p{e)de with p{e) being the renor- 
malized density of states, A the ultraviolet cutoff. For 
the spectrum we have the bare density of states 



Po(e) = (ttV* 2 ~ (2i') 2 ) 



hi- 



At 



■/■'■■I 



(4) 



which diverges logarithmically at e = —p. We will 
show that the renormalized density of states contains a 
much stronger divergence, which results in the pinning 
of the Fermi surface to VHS. To this end we calculate 
the electron self-energy S(k, e) for k near the VH point 
kvH — (tt, 0). In the second order in U this contains 
three contributions that arise from intermediate quasi- 
momenta which are close to the same or another VH 
point (cf. Refs. §,0), £( 2 )(k,e) = £- = i E *( k ' £ ) with 



Re£i(k,e) = - {2g / sin2p) 2 [A l£ + Bfk 
Re£ 2 , 3 (k£) 



Blkl]\n 2 {K/v) 



B x k 2 



2( 5o /sin2^) 2 [ J 4 2 , 3 £-f 
+2C 2 , 3 k x k y + Bl 3 k 2 y ] In(A/u) (5) 



where k x = 



e x sm ip, 



k y cos ip, g 



U/(4tt t) is the dimensionless coupling constant, 
m&yi(kl,k 2 y,\e\/2t,\p\/2t), A 1 = In 2, Bf = -Bf = 
1/2 — In 2. The analytical expressions for the coefficients 
Ai, Bi, and Cj with i — 2,3, which are some regular func- 
tions of t'/t, will be published elsewere. We perform a 
summation of logarithmically divergent terms which de- 
pend on A — (1/2) ln(A/u) within the RG approach. To 
separate effects of the electron spectrum renormalization 
and multiplicative renormalization of the Green's func- 
tion, we perform the renormalization procedure in two 



steps. At the first step we introduce the k-dependent 
mass renormalization factors Z^^v) (a = x,y), and at 
the second one all the other divergences are absorbed into 
the energy-dependent quasiparticle residue Z(e), 

G(k E ) = Z ( £ ) 

e + 2tZ x ^{v)kl - 2tZy^{v)kl + p + »T(e) 

(6) 

where Jl = p + Re£(ky# , —fi) is the renormalized chem- 
ical potential, the damping T(e) = Z (e)lmY,(kv h , e) is 
determined by analytical continuation. The coefficients 
Z am (v) = Z»j (v) Z^ 2 (v) Z^ 3 (v) satisfy the RG equa- 
tions 



dX 

d\nZ^ 2 (\) 
dX 

dAnZUX) 
dX 



= (B?-A 1 )Xjl/sm 2 2p 

= {B a 2 - A 2 )( 1 2 l + 72 2 - 7l72 )/ sin 2 2p 

= (5 3 a - A 3 ) 7 2 / sin 2 2p (7) 



Here 7, are four-electron vertices determined from 



dji/dX 
dj 2 /dX 
d^/dX 
dji/dX 

where 



4di(A)7i(7 2 - 71) + 4^7x74 - 4d 3 7!7 2 
2di(A)( 7 I + 7s 2 ) + 4d 2 (7i - 72)74 - 2d 3 ( 7l 2 + 7 |) 
-4d (A)7 3 74 + 4di(A)7 3 (27 2 - 71) 
-2d (A)( 73 2 + 74) + 2rf 2 ( 7 i 2 + 2 7 i72 ~ 2 7 2 + 74 2 ) (8) 



d (A)-A(l-i? 2 )- 1 / 2 ; 
di(X) = min(A,ln[(l + y/l -R 2 )/R}] 
d 2 = (l- R 2 y 1/2 ; d 3 = tan _1 (.R 



l-R 2 )/R (9) 



with R = -2t'/t. Equations for Z(v) = Z x {v)Z 2 {v)Z z {v) 
are obtained by the replacement Bf — Ai — » Ai in (0). Al- 
though the renormalization-group equations (j?]) and (j^) 
are only approximate because of the presence of ln 2 (A/u) 
terms in the series expansions, it was argued in Ref. |l3| ] 
that they reproduce correctly the solutions of the parquet 
equations which take into account complete momentum 
dependence of the vertices. Numerical solution of Eqs. 
(0) demonstrates that the vertices 7$ increase with de- 
creasing v and diverge at the critical value v = v c (t' /t) 
fill . We have Z am (v) — > 00 at v — ► v c while Z a (v) — > 0. 
This divergence signals a transition into magnetically or- 
dered or superconducting state at \p\ < p c — 2tv c . The 
equations (]?]), (|^) are valid in weak- and intermediate- 
coupling regimes where 7, < 1, i.e. at \p\ not too close 
to p c . 

The electron spectrum is determined by the pole of 
the Green's function (0). Increasing the factors Z am (v) 



—2 



with decreasing v leads to \s + fi\ < 2tmin(k x , k 2 ) at 



,-r2 



the pole, so that v = max(fc^, k 2 , \p\/2t) to logarith- 
mic accuracy. This demonstrates that the separation of 



2 



momentum- and energy-renormalization effects turns out 
to be self-consistent. We obtain therefore for the renor- 
malized spectrum 



have at <C \e\ <C t to leading logarithmic accuracy 
for £(ky#,e) the result 



2t 



k 2 

J) 



Z-^{\{i\/2t)k 2 y 



- Z-i(\ji\/2t)kZ 



k 2 

k 2 



k 2 CHAt/\e\) 



\jl\/(2tj r '~' n\e\ [l + Cln 2 (Ai/|e|)] 2 

\Jl\/(2t) 

Therefore we obtain 



C 



sin 2 2ip 



(10) 



For \n\ > /i c and k 3> k c 



,1/2 



we have 2 ' (* 2 ) 



and the dispersion coincides with the bare one, while at 
k x , k y < k c the spectrum becomes more flat due to renor- 
malization by the factors Z~^ lym (\Jl\/2t). Although the 
regime < fj, c , k < k c is outside of the region of the 
validity of RG equations ([j]) and (g), we have formally 
i?k = — M f° r k < k c , i.e. the spectrum is flat in this 
region. This is similar to the result of Dzyaloshinskii || . 
However, unlike that paper, we have taken into account 
all the channels of electron scattering, which gives a flat 
part of the spectrum already in one-loop approximation 
for the vertices. The flat part formation was considered 
earlier as one of formally possible instability channels of 
Landau Fermi- liquid ("fermionic condensation" 
Here we demonstrate that similar phenomena can take 
place in two-dimensional systems near the van Hove fill- 
ing. The result of numerical calculation of the dispersion 
law .Ek according to Eqs. (0) and (|ifj| ) at different values 
of Jl is shown in Fig.l. One can see that even at > fi c 
the spectrum has a wide almost flat region. The picture 
shown in Fig.l is in a good qualitative agreement with 
the ARPES data for La 2 _xSr x Cu0 4 These data 

give a possibility to estimate k c ~ 0.6a^ (a is the lat- 
tice constant) for this system; according to (|J), this k c 
value corresponds to ju c /(2i) = v c ~ 0.1. We stress once 
more the difference of the present approach with the ap- 
proaches [^3| which yield the flattening of the spectrum 
owing to strong antifferomagnetic correlations and there- 
fore are reliable only at low hole concentrations. 

The flat part of electron spectrum leads to drastic 
changes in the density of states and dependence n(JT). 
The contribution of the flat part of the spectrum to the 
density of states can be written as 



6p(e) = - 



k 2 . 



ImS(ky ff ,e) 



TT 2 [ e + Jl - ReS(k yff , e)] 2 + [Im£(k yH , e)] 2 

(11) 



k 2 

-^A(k VH ,s 

TT 



where £(k,e) = E(k, e) — ReS(k, — fi) and A(kvH,s) is 
the quasiparticle spectral weight at the van Hove momen- 
tum. 

Consider first the results obtained within the second- 
order expression for £(ky#,£), Eq. (||). Note that at 
this stage the effects connected with the renormaliza- 
tion of electron dispersion are already taken into account 
by performing summation of k-dependent logarithmically 
divergent terms to all orders of perturbation theory. We 



n(Jj) = n VH 



n 2 ^/W 



In — 

(2t>) 2 H 

sign Jl 



(12) 



(13) 



2Til + Clii 2 (At/\Jl\) 

where nyn is the van Hove filling. The second term in 
the right-hand side of Eq. ([13]) comes from momenta 
k > k c outside the flat part of the spectrum and can be 
calculated by using the bare density of states (||) since 
the renormalization of the spectrum at these momenta is 
not too important. For small enough \J1\ we can neglect 
this term in comparison with the third term to obtain 



n(n) = n VH 



k 2 sin 2ip sign \x 
2itg 2 \n2 \n 2 {kt/\]i\y 



(14) 



(the unity is small in comparison with squared logarithm 
in the denominator). The chemical potential measured 
from the VHS energy is given by a non-analytical func- 
tion, 



n(n) — At exp(— const/ \n — nval 1 ^ 2 ), 



(15) 



and is therefore practically zero in a rather wide range of 
electron concentrations n near uvh- Generally speaking, 
in the absence of the quadratic terms in the electronic 
spectrum, quartic terms may become important. How- 
ever, the presence of such terms give only subleading cor- 
rections to the above results and does not change them 
qualitatively. 

The results of numerical calculations of n{Ji) with Eq. 
( |l3| ) and the renormalized dependence n'{]l) with account 
of renormalization of £(ky#, e) according to Eq. (0) are 
shown in Fig. 2 for U — At and k c ~ 0.6. The depen- 
dence no(]i) obtained by integrating the bare density of 
states (Q) is also shown for comparison. The result for 
n'(jl) is shown at Jl < —fj, c only. One can see that the 
dependence n(jl) yields the pinning of the Fermi surface 
in the concentration range about 4% above and below 
VH filling (in fact, in the present approach the picture is 
symmetric, n(Jl) — tivh = n VH ~ n i~V))- Such a behav- 
ior is in a qualitative agreement with the ARPES results 
of Ref . [[l5| . On the other hand, this is in contrast with 
the behavior of chemical potential measured relative to 
its position in the insulating phase, Ref. j24|. 

It is important that the condition \n'(JT) — tivh\ > 
\n(JT) — uvh\ holds. This means that after account 
of renormalizations the pinning effect becomes even 
stronger than that given by Eqs.(p"3|) and ([ll]). Being 
extrapolated to the region \Jl/t\ < v c , the dependence 



3 



n'(Jx) should give even larger values of critical concentra- 
tions (about 6%). One can expect that for larger values 
oiU/t the critical concentrations will be larger. 

To conclude, wc have developed a RG approach which 
gives a possibility to describe the flattening of the elec- 
tron spectrum and does not suppose the presence of 
strong antiferromagnetic fluctuations. As a result of the 
peculiar structure of the spectrum, the pinning of the 
Fermi level to van Hove singularities in 2D systems oc- 
curs, the chemical potential being practically constant in 
a range of dopings near VH filling. In the pinning re- 
gion the electron density of states is determined by the 
quasiparticle damping and the system demonstrates es- 
sentially non-Fcrmi-liquid behavior. Further experimen- 
tal investigations on LaSrCuO and Bi2212 systems would 
be of interest to verify the pinning picture proposed. 
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Figure captions 

Fig.l. Quasiparticle dispersion for t'/t = —0.3 and 
U — At from RG approach. The values of the chemical 
potential are Ji = 0, — 0.2£, — OAt (from top to bottom). 

Fig. 2. The dependence n(Jl) — nyii for the same pa- 
rameter values as in Fig.l, Jl being referred to VHS en- 
ergy. The adopted value of k c is 0.6. The dashed line 
corresponds to bare electron spectrum, and the solid line 
to the result (|l^) . The dot-dashed line is the dependence 
with account of renormalizations of £(ky# , e). 
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